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Introduction: Several recent studies have pro- 

posed certain terrain types (i.e., physiographic units) 
on Titan thought to be formed by fluvial processes 
acting on local uplands of bedrock or in some cases 
sediment [e.g., 1-4]. We have earlier used our land- 
form evolution models to make general comparisons 
between Titan and other ice world landscapes (princi- 
pally those of the Galilean satellites) that we have 
modeled the action of fluvial processes [5]. Here we 
give examples of specific landscapes that, subsequent 
to modeled fluvial work acting on the surfaces, pro- 
duce landscapes which resemble mapped terrain types 
on Titan. 

Initial Terrain type Investigation: Crenulated 
Terrain and Isolated Ridges are commonly found in the 
equatorial latitudes of Titan, especially in the Xanadu 
region. Crenulated terrain generally occurs in discrete 
patches (sometimes closely spaced), often with signifi- 
cant linear elongation that might be associated with 
tectonic deformation [6]. Crenulated terrain has been 
interpreted by us and others to be mountainous relief 
that has been highly dissected by fluvial processes [1- 
6]. Some crenulated terrain is associated with dis- 
sected rims of large impact basins [3]. Such terrain 
may be analogous to eroded mountain belts on Earth. 
In the examples shown in Figure 1, crenulated uplands 
can be seen to range from nearly continuous in expo- 
sure to isolated belts surrounded and separated by in- 
creasing amounts of sediment infilling local lows. We 
investigate here whether the fluvial erosion of a region 
of tectonized terrain comprised of belts of ridge and 
trough (grooved) terrain separated from one another by 
regional lowlands might result in a landscape that re- 
sembles the crenulated and ridged terrains of Titan. Or 
alternatively would the erosion of a randomly- 
generated pseudo-fractal terrain of varying degree of 
small-scale roughness yield a closer match. 

Landform Evolution Modeling: We use a land- 
form evolution model (LEM) [7] to simulate fluvial 
and lacustrine modification to explore the degree to 
which strongly eroded landscapes retain information 
about their initial morphology. We utilized two types 
of initial conditions: (1) Digital elevation models 
(DEMs) of actual icy work landscapes not modified by 
fluvial processes (primarily tectonized landscapes of 
Ganymede and Europa) supplied by PMS; and (2) 
Randomly-generated pseudo-fractal terrain of varying 


degree of small-scale roughness. For type 1 runs, we 
assumed complete runoff, so that depressions become 
lakes, and drainage exits somewhere along edges of 
simulation domain. Edges of simulation domain were 
fixed (non-erodible). For type 2 runs, edge boundaries 
were periodic and runoff was balanced by lake evapo- 
ration, resulting in small lakes within the simulation 
domain. Flow rates were scaled to Titan gravity. Rate 
of channel incision was set to be proportional to shear 
stress on bed (parameterized as a power function of 
discharge and gradient). Weathering in the model pro- 
duced transportable debris. Sediment was deposited in 
low-lying areas, producing alluvial fans and infilling 
depressions from the edges inwards in deltaic deposits 
within lakes. We assumed little fine suspended sedi- 
ment deposited in lakes. In LEMs, the spatial pattern 
of erosion does not vary much with different parame- 
ters but the rate of erosion is very affected. For the 
purposes of these initial trials, the pattern not the rate, 
was deemed most important. 

Results: The pseudo-fractal landscapes have 

roughness expressed over multiple spatial scales, but 
no expression of systematic structural influence (Fig. 
2a). When eroded, highlands are eroded by weathering 
followed by fluvial entrainment, whereas lowlands are 
infilled by bajada-like fans and deltaic deposits gran- 
dually encroaching on lacustrine depressions (Fig 2b). 
Because of the balance between lacustrine evaporation 
and runoff, persistent lakes are formed which deepen 
through time, but they may become inactive and in- 
filled if drainage is subsequently integrated into larger 
basins. If deposition of suspended sediment were mod- 
eled, these depressions would feature infilling by fine 
lacustrine sediment. The landscape becomes bimodally 
segregated into steep, dissected uplands with abrupt 
but irregular lateral transitions into smooth, low gradi- 
ent alluvial plains (Fig. 2b). The location of highlands 
and plains is dictated by the original topography. If the 
landscape is subjected to further erosion, the uplands 
shrink and lower, and the alluvial plains expand to 
cover a greater percentage of the landscape. 

Simulated erosion of tectonized terrain from Ga- 
nymede (30°S, 185°W, Fig 3a) similarly results in 
three landform elements, steep eroded uplands, broad 
alluvial plains, and depressions gradually infilled by 
deltaic deposits, particularly evident in the interior of 
craters, which, in this simulation, are assumed to host 



deep lakes due to the complete runoff and no assumed 
evaporation. Particularly noteworthy is the strong in- 
heritance of initially elevated regions in the resulting 
eroded landscapes, which, in this case, features linear 
ridges inherited from the grooved initial surface (Fig. 
3b). 

Discussion and Conclusions: (1) Regarding the 
possible structural control of eroded ridges on Titan. A 
surprising result from this study is that eroded pseudo- 
fractal terrain produces a landscape broadly composed 
of belts of isolated somewhat aligned ridges. As noted 
above, the erosion of the Ganymede sample produces a 
landscape in with the individual ridges exhibiting a 
higher degree of alignment and similarity of scale than 
is observed on Titan. We are not finding eroded ex- 
amples of Ganymede- style “grooved terrain” on Titan. 
(2) Alluvial lowlands in fractal terrain, if partially 
flooded by increase in methane runoff, would resemble 
lakes on Titan - if shallowly flooded, would be like 
Ontario Lacus - if more deeply flooded, like the N 
polar lakes with some portions of shorelines transition- 
ing to alluvial plains and others directly abutting 


eroded highlands. Lakes could shallow, as observed 
[8], despite the abrupt, crenulated boundaries with 
eroded uplands. (3) If erosion is continued until few 
uplands remain, our models appear to produce possible 
analogs to dissected mottled terrain (allowing for some 
of this terrain to feature aeolian mantling as well). We 
are currently beginning runs to evaluate the original 
appearances of dissected mottled terrain, putative 
cratered terrain, and dissected plateau terrain [1]. 
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Fig. 1. Three views of equatorial Crenulated Terrain or Isolated Ridges [1] interpreted to be exposures of uplands 
bedrock, principally eroded by fluvial activity, showing increasing degrees of local sediment infilling of the local 
lows (left to right). Individual frames cover an area -275 km across, N up. All at ~10°S. a. 1 15°W, b. 100°W, c. 
205°W. (Excerpts from the latest USGS digital base map.) 



Fig. 2. The pseudo-fractal land- 
scapes have roughness expressed 
over multiple spatial scales, but no 
expression of systematic structural 
influence. A surprising result from 
this study this that fluvial erosion 
of this terrain produces a landscape 
broadly composed of belts of iso- 
lated somewhat aligned ridges. 



Fig. 3. Simulated fluvial erosion of 
DEM of Ganymede terrain (-300 
km across). The erosion of the Ga- 
nymede sample produces a land- 
scape in with the individual ridges 
exhibiting a higher degree of 
alignment and similarity of scale 
than is observed on Titan. 



